Effect of nonmagnetic disorder on criticality in the "dirty" U(l) spin liquid 
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We investigate the effect of nonmagnetic disorder on the stability of the algebraic spin liquid (ASL) 
by deriving an effective field theory, nonlinear a model (NLaM). We find that the anomalous 
critical exponent characterizing the criticality of the ASL causes an anomalous gradient in the 
NLaM. We show that the sign of the anomalous gradient exponent or the critical exponent of 
the ASL determines the stability of the "dirty" ASL. A positive exponent results in an unstable 
fixed point separating delocalized and localized phases, which is consistent with our previous study 
[Phys. Rev. B 70, f 40405 (2004)]. We find power law suppression for the density of spinon states 
in contrast to the logarithmic correction in the free Dirac theory. On the other hand, a negative 
exponent destabilizes the ASL, causing the Anderson localization. We discuss the implication of 
our study in the pseudogap phase of high T c cuprates. 

PACS numbers: 71.30.+h, 72.15.Rn, 71.23.-k, ll.10.Kk 



03 



I 

C 

O 
o 



(N 
> 

in 

o 
in 
o 

03 



o 
o 



X 



The effect of nonmagnetic disorder in the presence of 
strong correlations between electrons is one of the central 
interests in modern condensed matter physics. This situ- 
ation is realized when nonmagnetic impurities are doped 
in the Mott insulator. For example, Zn doped high T c 
cuprates are those. In order to solve this problem it 
should be clarified how to describe the Mott insulator. 
Here we concentrate on a paramagnetic Mott insulator 
(PMMI). According to one possible scenario, hole dop- 
ing to a PMMI results in high T c superconductivity 
Understanding the nature of the PMMI may be cru- 
cial for the mechanism of high T c superconductivity. Re- 
cently, the PMMI was proposed to be the U(l) spin 
liquid (t/lSX)0. The U1SL is the state described by 
QED3 in terms of massless Dirac spinons strongly in- 
teracting via noncompact U(l) gauge fields 2]. The key 
feature of the U1SL lies in the criticality characterized 
by anomalous critical exponents originating from long 
range gauge interactions [1 H H H H Ufl- All cor- 
relation functions show algebraic power law decay with 
anomalous critical exponents,resulting in no well defined 
quasiparticles (spinons) [1 H H If! H H IqL ITci| . Owing to 
the power law behavior the U ISL is sometimes called the 
algebraic spin liquid (ASL) more appropriately. 

Recently, the role of nonmagnetic impurities in 
the ASL was investigated by the present author [TT| . 
In contrast to the free Dirac theory in two spacial 
dimensions|l2l Il3| long range gauge interactions can in- 
duce a delocalized state at zero temperature, .llj. The 
presence of nonmagnetic disorder destabilizes the free 
Dirac fixed point. The renormalization group (RG) flow 
goes away from the fixed point, indicating localization|12L 
. On the other hand, the ASL fixed point in QED^M 
remains stable at least against weak disorder llj. An 
unstable fixed point separating delocalized and localized 
phases is found The RG flow shows that the effect 
of disorder potential vanishes if we start from sufficiently 
weak disorder. 

In the present study we reexamine the effect of non- 
magnetic disorder on the ASL by deriving an effective 



field theory, nonlinear a model (NLaM) . The main idea 
is firstly to integrate out gauge fluctuations. This gives 
the single particle propagator Eq. (2) with the anoma- 
lous exponent 77^,, the hall mark of the ASL. Based on 
this ASL propagator we obtain the NLaM Eq. (8) by 
usual treatment. The core in the NLaM is the anoma- 
lous gradient exponent rj a originating from the critical 
exponent 77^ of the ASL. We find that the sign of r\ a 
or riif) determines the stability of the ASL against non- 
magnetic impurities. In the case of r\ a > an unstable 
fixed point separating delocalized and localized phases is 
found, which is consistent with our previous study [Tl|. 
The density of spinon states and the conductance of a 
spinon number (internal gauge charge) show power law 
corrections in contrast to the logarithmic suppression in 
the free Dirac theory [13|. In the case of 77^ < the 
ASL becomes unstable. The Anderson localization is 
expected. We discuss the implication of our study in the 
pseudogap phase of high T c cuprates. 



First we briefly review how the effective QED3 La- 
grangian is derived from the antiferromagnetic Heisen- 
berg model, H — JJ2<ij> &i ' &j wu; h J > 0. Insert- 
ing the fermionic representation of spin Si = \ f} a T a p \fi@ 
into the Heisenberg model, and performing the stan- 
dard Hubbard-Stratonovich transformation for an ex- 
change hopping channel, we obtain an effective one body 
Hamiltonian for the fermions coupled to an order pa- 
rameter, H ef f = -JE<i,i> fLxijfjc* ~ h.c. Here 
fia is a fermionic spinon with spin a =|, j, and Xij 
is an auxiliary field called a hopping order parame- 
ter. Notice that the hopping order parameter Xij is a 
complex number defined on links ij. Thus, it can be 
decomposed into Xij = \Xij\ el6ij where \xij\ and #ij 
are the amplitude and phase of the hopping order pa- 
rameter, respectively. Inserting this representation of 
Xij into the effective Hamiltonian, we obtain H e ff = 
~ J Y«i,j > \Xij\fi a el0ii fj a -h.c. Then, we can easily see 
that this effective Hamiltonian has internal U(l) gauge 
symmetry HLJf^dU = H eff [f ia ,e ij ) under the fol- 
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lowing U(l) phase transformations, f' ia — e i ^ >i fi a and 
Q'ij = % + & — fy- This implies that the phase 9^ of 
the hopping order parameter plays the same role as the 
U(l) gauge field Oy. When a spinon hops on lattices, 
it obtains the Aharnov-Bohm phase owing to the U(l) 
gauge field ay. It is well known that the stable mean 
field phase is ir flux state at half filling[f|. This means 
that a spinon gains the phase of n when it turns around 
one plaquette. In the n flux phase low energy elemen- 
tary excitations are massless Dirac spinons near nodal 
points showing gaplcss Dirac spectrum and U(l) gauge 
fluctuations^!- In the low energy limit the amplitude 

\Xij\ is frozen to be \xij\ = J\ < /j Q /«* > I = Xo- As 
a result we obtain the following low energy effective La- 
grangian in terms of massless Dirac fermions near the 
nodal points interacting via compact U(l) gauge fieldsQ 
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Here tp ncr is the two component massless Dirac spinon, 
where n = 1,2 represent the nodal points, (ir/2,ir/2) 
and (tt/2, -tt/2), and a =|, |, SU(2) spin. They are 

given by ^i<r = ( { lecr ) and ^2<r = ( { 2ocr ) , respec- 

\ hoa J \ J2ea J 

tively. In the spinon field f n i& n — 1,2 represent the 
nodal points, I = e,o, even and odd sites, and a =|, J,, 
its spin, respectively p|- The Dirac matrices 7^ are given 
by the Pauli matrices 7^ = (0-3, 02, o"i) where they satisfy 
the Clifford algebra [7^,7^] + = 2<5 /i „[j|. is the U(l) 
gauge field. The kinetic energy of the gauge field results 
from particle-hole excitations of high energy spinons. e 
is an effective internal charge, not a real electric charge. 

We have difficulty in solving Eq. (1) owing to instan- 
ton excitations originating from the compactness of the 
U(l) gauge field. Instantons represent tunnelling events 
between energetically (nearly) degenerate but topologi- 
cal^ incquivalent vacua. In the U(l) gauge theory the 
instanton has a magnetic monopole configuration. Mag- 
netic monopoles are known to play a crucial role in con- 
finement physics^, El- Recently, it was shown that the 
instanton effect can be suppressed at least in the large 
N limit (a =|, 4 o = 1, 2, TV) where N is the flavor 
number of massless Dirac fermions 0. Existence of quan- 
tum criticality in the QED3 with noncompact U(l) gauge 
fields is the key mechanism of suppression of instanton 
excitations. It is well known that the QED3 has a stable 
charged fixed point in the large N limit. Hermele et. al 
examined the stability of the charged fixed point against 
instanton excitations 2]. At the fixed point charges can 
be sufficiently screened by critical fluctuations of Dirac 
fermions. The larger the flavors N of the Dirac fermions, 
the smaller electric charges, thus resulting in huge mag- 
netic charges. This can cause a negative scaling dimen- 
sion to instanton fugacityQ. As a result instanton ex- 
citations can be suppressed at the charged fixed point 



in the large N limit. Deconfinement of charged matter 
fields, here the massless Dirac spinons can be achieved 
at the quantum critical point. The QED3 with noncom- 
pact gauge fields appears as the critical field theory at 
the charged fixed point. All correlation functions exhibit 
power law behaviors with anomalous critical exponents 
generated from long range gauge interactions. In this re- 
spect the state described by the QED3 is called the ASL 
and the charged fixed point, the ASL critical point. 

As discussed in the above, the quantum criticality of 
the effective QED 3 Eq. (1) is the crucial point for de- 
confinement of the Dirac spinons. But, we should re- 
member that the ASL criticality can survive only in the 
limit of large flavors N corresponding to the SU(N) quan- 
tum antiferromagnet (a — 1,2,..., TV). In the case of 
the physical SU(2) antiferromagnet (a it is not 

clear if the ASL criticality remains owing to sponta- 
neous chiral symmetry breaking (SxSB) causing anti- 
ferromagnetism. It is believed that there exists the crit- 
ical flavor number 7V C associated with the SxSB in the 
g^^HESEm. But > the precise value of the criti- 
cal number is far from consensus|l6j|. If the critical value 
is larger than 2, the SxSB is expected to occur for the 
physical N = 2 case. Then, the Dirac fermions become 
massive. In this case the massive Dirac spinons are con- 
fined to form spin 1 excitations, antifcrromagnons0. On 
the other hand, in the case of N c < 2 the ASL criticality 
remains stable against the SxSB. In the present paper 
we assume that the chiral symmetry is preserved even in 
the SU(2) quantum antiferromagnet. There is a support- 
ing argument for existence of the ASL even in the case 
of N c > 20. If there are additional massless fluctuations 
carrying internal U(l) gauge charges, the SxSB can be 
forbidden. These additional gapless fluctuations can arise 
from quantum critical points in association with antifer- 
romagnetism or superconductivity. For example, criti- 
cal fluctuations of doped holes near the superconducting 
transition quantum critical point are expected to increase 
the flavor number of massless fluctuations |12|. If the to- 
tal flavor number of massless Dirac spinons and doped 
holes exceeds the critical value N c , the SxSB is not ex- 
pected to occur [Til ITsf . In this case the quantum critical 
point would be described by the ASL for spin degrees of 
freedom lj. The role of doped holes in the ASL will be 
discussed later in more detail. 

The criticality of the ASL is characterized by the crit- 
ical exponents of correlation functions. The single par- 
ticle propagator GASh{k) =< Vw(&)Vw(fc) > is given 
byHHHHEa 



G AS L{k) 



■ 7 A' 



(2) 



where 77^, is the anomalous critical exponent of the 
ASL propagator. Here we briefly sketch how the 
ASL propagator Eq. (2) is derived from the effec- 
tive QED 3 Eq. (1). The single particle propaga- 
tor is generally given by G A l SL {k) = G' G 1 {k) + E(fc), 
where Gq (k) = ij^k^ is the inverse of the bare spinon 
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propagator, and the spinon self-energy result- 

ing from long range gauge interactions. In the usual 
1/N expansion^, IE 01 the self-energy is represented 
by S(fc) = / jffisl^Goik + qhvD^iq), where D^(q) 

is the renormalized propagator of the U(l) gauge field 
due to particle-hole excitations (polarization) of mass- 
less Dirac fermions. The gauge propagator is obtained to 

in the Lorentz 



be D^(q) «n- x (g) = ^(& 



gags 
q' 2 



gauge, where Up V (q) 



is the polarization function of the Dirac fermions. In- 
serting this gauge propagator into the expression of the 
self-energy, we obtain the spinon self-energy of loga- 
rithmic momentum dependence E(fe) = irj^^k^lnf^j , 

where is the anomalous exponent and A, the momen- 
tum cutoff. The absolute value of the exponent xty is 
proportional to the inverse of the flavor number, i.e., 

Iwi ~ n -1 0,111111 OH Since the QED3 is the critical 

field theory at the charged critical point, all correlators 
should exhibit power law behaviors. In this respect the 
logarithmic momentum dependence should be considered 
to be the lowest order in the power law behavior. As a re- 
sult we can obtain the ASL single particle propagator E< 
(2) from the following nonperturbative consideration 



G ASL( k ) 
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though the spinon propagator exhibits its algebraic form, 
it is difficult to give a definite physical meaning. This is 
because it is not gauge invariant. All physical observ- 
ablcs should be gauge invariant. The critical exponent 
770 should be evaluated in a gauge invariant way. The 
following gauge invariant green function is usually con- 
sidered, Gasl(x) =< T T [^ na {x)e^ ^^(0^(0)] >. 
It is not easy to obtain the critical exponent 770 by cal- 
culating this gauge invariant green function. Its precise 
value is far from consensus and under current debate. 
The crucial point in the following discussion is the sign 
of the exponent rjy,. Most evaluations [5. Isl IsL ITo| suggest 
its negative sign, < 0. However, as argued in Ref. [6|, 
its negative sign is unphysical in the sense that the ASL 
propagator is more coherent at long distances than the 
propagator of the free Dirac theory owing to long range 
gauge interactions. This result is completely in contrast 
to the usual role of interactions. Interactions would make 
the propagator less coherent. This is indeed true in the 
critical field theories with local repulsive interactions (for 
example, <^> 4 theory). Positive critical exponents are well 
known in these theories. If the critical exponent r/^p is 
positive, the long range gauge interactions destabilize the 
Fermi liquid pole (the pole of the single particle green 
function in the free Dirac theory). The renormalization 
factor Z(p) representing weight of quasiparticles is given 
by Z(p) ~ p 71 ^ with momentum p. In the case of f]^ > 
it vanishes in the long wave length and low energy limit, 
i.e., Z(p — > 0) — > 0, leading to Luttinger liquid-like power 
law correlators. This is the ASL as a critical state. In 
this respect the ASL can be considered to be two di- 
mensional realization of the one dimensional Luttinger 



liquid|(|. In the present paper we do not determine its 
sign. Instead we use the exponent 770, as a phenomeno- 
logical parameter. We consider both cases, 770 < and 
770 > 0. Here we assume that the ASL green function 
Eq. (2) is obtained in a gauge invariant wav|3.l^.l8ll9lllo| 
and thus the critical exponent rj^, is gauge invariant. We 
repeat that its absolute value isgiven by \rj^\ ~ A^" 1 in 
the 1/N approximation |3l 161 ISL l9l Il0| . 

The ASL propagator G A SL(k) =< Vw(fc)Vw(fc) >= 
—i fc 1i-t^ can De easily obtained from the following La- 
grangian 



(3) 



Integration over U(l) gauge fluctuations in the QED3 
Eq. (1) results in the ASL green function Eq. (2) as 
explicitly demonstrated above. Eq. (2) is easily derived 
from Eq. (3). Therefore, Eq. (3) can be considered to 
be an effective field theory resulting from integration over 
the U(l) gauge field in the QED3 Eq. (1). In this respect 
the Dirac spinon field ip nr7 in Eq. (3) is clearly different 
from that in Eq. (1). It should be considered to be a 
renormalized field resulting from the gauge interactions. 
The renormalized Dirac spinon field arises from the self- 
energy correction. We view Eq. (3) as the prototype 
describing the ASL. More generally, one may think that 
Eq. (3) represents one class of critical field theories de- 
pending on the critical exponent 770 . This point of view 
is parallel to the standpoint that a free fermion theory 
is the foundation describing Fermi liquid. One cautious 
theorist can argue that the ASL Lagrangian Eq. (3) is 
not sufficient because Eq. (3) does not include appro- 
priate vertex corrections. But this guess is not correct. 
Remember that the critical exponent 770 was obtained in 
a gauge invariant way although the calculation of 770 in 
a gauge invariant way was not explicitly demonstrated 
in the above owing to its complexity pi. la. Isl Wi ITfl| . It is 
well known that we cannot satisfy the gauge invariance 
without vertex corrections 0- Thus, the evaluation of 
770 in a gauge invariant way includes vertex corrections. 
This argument can be checked by considering a two par- 
ticle green function. The critical exponent of a spin-spin 
correlation function, obtained in a gauge invariant way 
including vertex corrections explicitly in Eq. (1), is given 
by twice the exponent of the single particle propagator 
in Eq. (3), i.e., 2rj^ in the case of 770 < 0|J,l3. This was 
also pointed out in Ref. [l(j- This implies that we can 
treat the ASL effective Lagrangian Eq. (3) as a "free" 
theory. 

We introduce a random potential V(r) in the ASL [Eq. 
(3)] 



S = 



J d ° x £ £ [^wrM^^Vw + V{v)i> n 



As the randomness is independent of time and the Casl 
is quadratic in ip na , each frequency sector in Eq. (4) is 
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separated from each other. Thus it is sufficient to con- 
sider only zero frequency sector. This motivates us to 
explore D = 2. Here we assume that V(r) is a gaussian 
random potential of < V(r)V(r') >= W5(r — r') with 
< V(y) >= 0. Using the standard replica trick to av- 
erage over the gaussian random potential, we obtain the 
following effective action in the presence of nonmagnetic 
disorder 

„ M 

S = j d 2 r^2 (* a T i a i 8-'»* a + ie^ a y c 



w M 



a,a f — l 



(5) 



Here * Q 




is the eight component spinor 



where a = 1 , M is an replica index and the limit 
M — > is to be taken at the end. = I ® er 3 ® 7^ 
is the eight-by-eight gamma matrix. / acts on SU(2) 
spin space and a 3 , different nodal points. We have in- 
cluded an infinitesimal imaginary potential e in order to 
generate correlation functions [13j. In Eq. (5) one can 
easily read the bare scaling dimension of the disorder 
strength W . It is given by dim[W] = —2rj^. In the 
case of m > this is consistent with our previous RG 
study [III]. In the study the RG equation is found to 
be dW/dlnL w —x^W to the first order in W. \ is 
a positive numerical constant. Considering the charged 
fixed point e 2 ~ N~ l with the flavor N, the RG equa- 
tion yields dimflU] = — x e c ~ — -W -1 ~ — ?/ty- But, in 
the opposite case of 77^, < the scaling dimension of the 
disorder strength W is positive, indicating instability of 
W = 0. 

Performing the standard Hubbard-Stratonovich trans- 
formation and the gaussian integration for the spinon 
field we obtain the following effective action in terms 
of the order parameter field Q 



'e// 



d 2 r 



1 

2W 



Tr[Q 2 (r)] 



Trln(T t d t d~^' + iQ 



(6) 



Here Q is the 8M x 8M matrix field. Its saddle point 
is given by Q aa ' — W < ^ a ^a' >■ The saddle point 
equation is obtained to be 



1 = 8W 



d 2 k 



1 



(2tt) 2 |k| 2 - 2 '^ + Q 2 



WF(Q 2 ). (7) 



Here we replaced Q QQ < with Q QQ < = S aa >Q. Since F(Q 2 ) 
is a monotonically decreasing function from infinity to 
zero in the case of 77,/, < 1, there exists only one nonzero 
solution for The relation of \r)^\ ~ TV -1 guarantees 

this consideration. In the free Dirac theory (77^ = 0) 
Q = Ae _7r / 4W/ with the momentum cutoff A is obtained. 



In the ASL the meaning of nonzero Q is not clear. In 
Fermi liquid it corresponds to the imaginary part of the 
self energy and thus the finite value of Q generates a finite 
density of quasiparticle states due to disorder. But, in 
the present case there are no well defined quasiparticlcs. 
Thus we cannot define the density of quasiparticle states. 
Mathematically, it gives an imaginary part to the single 
particle propagator. In this respect we can think that a 
diffusive behavior of critical spinons emerges by disorder. 

Next we consider small fluctuations around this saddle 
point. Here we follow Ref. [2(j. Inserting the expression 
Qaa' = QU aQ / into Eq. (6) and expanding the loga- 
rithmic term to the second order in \5 aa ' , we obtain the 
NLaM^ 

Snl*m = J d 2 r(^-Tr[VUt(r)^-VU(r)] 
+eTr[U(r)+U t (r)]). (8) 

Here g~ l is considered to be the conductance a sn of a 
spinon number (internal gauge charge). Its bare value 

g%~ is obtained to be g° a ~ = | / ^[ fcM -*% +Q3 • 
The anomalous gradient exponent rj a is given by ?7 CT = 
277^,. As mentioned in the introduction, the anomalous 
gradient is the crucial feature in the NLaM . Its exis- 
tence can be interpreted to be the mirror of the critical- 
ity 77^, in the original system, here the ASL. Generally 
speaking, we are considering the role of nonmagnetic dis- 
order in the critical systems with anomalous critical ex- 
ponents. The stability of the ASL critical point is de- 
termined by the anomalous gradient exponent 770- or the 
critical exponent 77^ of the ASL. In the case of 77^ > 
(77^, > 0) the ASL critical point remains stable against 
weak randomness. But in the case of 77^ < (77,/, < 0) 
it becomes unstable and the Anderson localization oc- 
curs. We discuss this main issue by investigating the RG 
equation of the stiffness parameter g a . 

One can easily derive the following RG equation in one 
loop order 



dg a 
dlnL 



= -(d-2 + Va )g <7 + Agl. 



(9) 



Here d is the spacial dimension and d — 2 is the present 
case. A is a positive numerical constant. In the case of 
77^, — > the above RG equation is reduced to that of the 
free Dirac theory Furthermore, this equation corre- 
sponds to Eq. (7) in our previous study 11]. Inserting 
g a — > W and r\ ~ N -1 — > e 2 into Eq. (9), we obtain 
the following RG equation of the disorder strength W, 
dW /dlnL = —e 2 W + AW 2 . This completely coincides 
with Eq. (7) in the previous study [Til]. 

Now we discuss the phases from the RG equation (9). 
First we consider the case of a positive exponent, 77^ > 0. 
Then, this RG equation shows an unstable fixed point 
g* = rid A consistent with our previous study [Tl|. In the 
case of weak randomness g a < g* the g a goes to zero 
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owing to the first term, implying the emergence of a de- 
localized spinon state. In the case of strong randomness 
g a > g* the g a goes to infinity owing to the second term, 
indicating localization of the spinons. Solving the RG 
equation (9), we obtain the spinon conductance 



(10) 



g® is the bare spinon number conductance 



where <r° 

and l e , the elastic mean free path|13|. The key feature 
is the power law correction resulting from the anoma- 
lous critical exponent rj^. This is in contrast to the free 
Dirac theory. In the limit r/ a — > the present conduc- 
tance formula is reduced to the logarithmic suppression, 

a sn = a® n — Aln^Y^j consistent with that of the free 

Dirac theory 01 . The density of the spinon states can be 
easily calculated in the delocalized regime. Although it is 
not clearly defined in the ASL, the terminology is used 
below as it is. The formal expression of the density of 
states is given by p = lim M ^o jffj < Tr[W + U] >[T^. 
We find that the density of states also has the power law 
correction 



P- Po 
Po 



1 

1" 



(l-^-L-^), 



(11) 



where B is a positive numerical constant. This expression 
shows a finite density of spinon states in the limit of L — > 
oo, given by p = po(l — Bg cr rj~ 1 l~ ri <'). In the limit rj a — > 
we also recover the logarithmic suppression p ~^° = 

—Bg^lnyj-^ in the free Dirac theory, 1.3] . On the other- 
hand, in the case of a negative exponent r\ a < the RG 
equation (9) shows a runaway characteristic for the g a . 
This implies the Anderson localization of spinons and 
thus the ASL disappears in the presence of nonmagnetic 
impurities. 

Next we discuss how the delocalized state of spinons 
can emerge. In two spacial dimensions the NLaM based 
on the free Dirac theory does not have a stable fixed 
point^l indicating a delocalized state. But, the pres- 
ence of a topological term such as a Wess-Zumino-Witten 
(WZW) term or Berry phase (9) term can result in a 
stable critical point |2l|. This is well known in the disor- 
dered metal and antifcrromagnetic spin chain with spin 
1/2. On the other hand, the derealization in the present 
study has nothing to do with such topological terms. 
It originates from the criticality of the ASL. It is well 
known that the NLaM has an unstable fixed point above 
two dimensions The criticality of rj^, > leads the 
ASL to be in (2 + 2rj^,) dimensions effectively. As a result 
the derealization emerges against weak randomness. We 
would like to stress that the mechanism of derealization 
in the critical phase (ASL) totally differs from that aris- 
ing from the topological terms. This is reflected in the 
RG equation. In the derealization induced by the topo- 
logical terms the stable fixed point appears from the g^ 
term in the RG equation [2l|. This is associated with de- 
structive interference effect of interactions. On the other 



hand, the derealization driven by the criticality arises 
from the g a term [Eq. (9)]. This is due to increase of 
effective dimensionality, as discussed above. 

If the WZW term appears in the NLaM Eq. (8), a 
stable strong coupling fixed point is expected to exist |2l| 
in both cases, 77^ > and 77^ < 0. The density of states 
would exhibit a power law behavior near the fixed point. 
Moreover, it is expected to vanish as energy decreases 
down to zero (2]]]. This is in contrast to the present case. 
However, it is not clear whether this new fixed point is 
stable against instanton excitations. The spinons would 
be renormalized by the disorder effect. It is difficult to de- 
termine how the renormalized spinons affect the internal 
gauge charge. This problem will be an important future 
work on the stability of the ASL against both instantons 
and disorder. 

Our present study has important implications in the 
role of nonmagnetic disorder in high T c cuprates. Ac- 
cording to one scenario^, the pseudogap (PG) state is 
proposed to be the ASL. Then, the ASL should be sta- 
ble against disorder because all samples include disorder. 
The stability of the ASL against disorder depends on the 
sign of the critical exponent rj^. As pointed out earlier, 
most evaluations [III II US support a negative critical 
exponent. This would destabilize the ASL against dis- 
order, resulting in the Anderson localization. A positive 
critical exponent 6] is necessary for the stability of the 
ASL. Our investigation requires more careful determina- 
tion of the critical exponent rj^. Furthermore, the present 
study tells us that the previous studies of disorder effects 
in the PG state l22l a re difficult to be applied. This is be- 
cause the studies |22j| are based on the free fermion theory 
ignoring strong gauge fluctuations. The existence of the 
delocalized state in the impurity doped ASL seems to be 
inconsistent with experiments 22] . In experiments a non- 
magnetic impurity is believed to localize spin 1/2. This 
trapped spin acts as a free spin, showing the Curie- Weiss 
behavior in the spin susceptibility p2l| . We expect that 
this inconsistency may be resolved by considering strong 
disorder. We note that strong disorder causes the local- 
ization of spinons. Zn impurity can be considered to be 
a strong scatterer owing to its compact electronic shell 
structure. 

In the present study we did not consider the effect of 
hole doping. Doped holes are represented by holons in 
the context of the slave boson theory. Holons can af- 
fect gauge fluctuations. For example, when holons are 
condensed, gauge fluctuations become massive via the 
Anderson-Higgs mechanism. As a result a free Dirac 
theory is obtained. The spinon can be localized even in 
weak disorder. The Curie- Weiss behavior can be easily 
understood based on the free Dirac theory in the super- 
conducting state [23|. In the PG phase the holons are not 
condensed. Gauge fluctuations remain massless. In this 
case damping effect in gauge fluctuations is expected to 
arise from holon contributions |24l |25| . The role of disor- 
der in the presence of dampedgauge fluctuations would 
be an interesting future work|25|. 
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Lastly, we would like to comment about the present 
approach based on the NLaM . The NLaM approach 
has some advantages, compared with the previous access 
based on the fcrmionic action llj. First, the NLaM ap- 
proach is more efficient than the fcrmionic one in investi- 
gating physics of phase transitions. The NLaM assumes 
the presence of a local order parameter. On the other 
hand, the previous study[lll is based on the absence of 
local ordering. If we start from a disordered phase with- 
out local ordering and approach a critical point associ- 
ated with a phase transition, a local order parameter is to 
emerge, strongly fluctuating. In order to describe critical 
fluctuations of the order parameter, we explicitly intro- 
duce the local order parameter and obtain an effective 
field theory of the order parameter field by integrating 
over the fermions. The resulting effective field theory is 
the Ginzburg-Landau free energy. As a matter of fact 
the NLaM is nothing but the Ginzburg-Landau free en- 
ergy formulation. If we start from the fermion action, 
we realize that it is not easy to describe phase transi- 
tions. This is because phase transitions are basically 
nonperturbative phenomena not captured by perturba- 



tive calculations based on the fermionic action. In order 
to describe the phase transitions it is necessary to sum 
infinitely many diagrams in the fermionic action. On 
the other hand, in the Ginzburg-Landau formulation the 
phase transitions are easily described. This is the rea- 
son why the NLaM is called an effective field theory 
for phase transitions. Second, it is also easy to calcu- 
late some physical quantities like conductance and den- 
sity of states. These quantities directly appear in the 
NLaM . We can easily obtain the conductance and den- 
sity of states as a function of energy and size of a system 
by the present RG calculation. Last, topologically non- 
trivial excitations in the NLaM are not captured in the 
fermionic action although these are not intensively dis- 
cussed in the present paper. 

In summary, we examined the effect of nonmagnetic 
disorder on the ASL criticality characterized by its 
anomalous critical exponent. In the case of a positive ex- 
ponent the critical point remains stable at least against 
weak disorder. But, in the opposite case the critical point 
becomes unstable. 
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